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I. INTRODUCTION 



The matter created in relativistic heavy-ion colhsions manifests a strongly collective hy- 
drodynamic behaviour ^ which is particularly evident in studies of the so-called elliptic flow 
[3]. A hydrodynamic description requires, strictly speaking, a local thermal equilibrium and 
experimental data on the particle spectra and elliptic flow suggest, when analysed within 
the hydrodynamic model, that an equilibration time of the parton^ system produced at the 
collision early stage is as short as 0.6 fm/c j3|. Such a fast equilibration can be explained 
assuming that the quark-gluon plasma is strongly coupled However, it is not excluded 
that due to the high-energy density at the early stage of the collision, when the elliptic flow 
is generated |i5|, the plasma is weakly coupled because of asymptotic freedom. Thus, the 
question arises whether the weakly interacting plasma can be equilibrated within 1 fm/c. 

Models that assume that parton-parton collisions are responsible for the thermalization 
of weakly coupled plasma provide a significantly longer equilibration time. The calculations 
performed within the 'bottom-up' thermalization scenario where the binary and 2 3 
processes are taken into account, give an equilibration time of at least 2.6 fm/c To 
thermalize the system one needs either a few hard collisions of momentum transfer of order 
of the characteristic parton momentum^, which is denoted here as T (as the temperature 
of equilibrium system), or many collisions of smaller transfer. As discussed in e.g. 0, the 
inverse equilibration time is of order gHn{l/g) T (with g being the QCD coupling constant) 
when the binary collisions are responsible for the system's thermalization. However, the 
equilibration is speeded up by instabilities generated in an anisotropic quark-gluon plasma 
jy^lIOl; as growth of the unstable modes is associated with the system's isotropization. The 
characteristic inverse time of instability development is roughly of order gT for a sufficiently 
anisotropic momentum distribution 1^, 11,0,0, 3- Thus, the instabilities are much 
'faster' than the collisions in the weak coupling regime. Recent numerical simulation (1^ 
shows that the instabilities driven isotropization is indeed very efficient. 

The isotropization should be clearly distinguished from the equilibration. The instabil- 
ities driven isotropization is a mean-field reversible phenomenon which is not accompanied 
with entropy production 0|. Therefore, the collisions, which are responsible for the 
dissipation, are needed to reach the equilibrium state of maximal entropy. The instabilities 
contribute to the equilibration indirectly, shaping the parton momenta distribution. And re- 
cently it has been argued 3| that the hydrodynamic collective behaviour does not actually 
require local thermodynamic equilibrium but a merely isotropic momentum distribution of 
liquid components. Thus, the above mentioned estimate of 0.6 fm/c jsj rather applies to the 
isotropization than to the equilibration. 

My aim here is to review the whole scenario of instabilities driven isotropization and 
the article is organized as follows. I start with a brief presentation of numerous efforts to 
understand the equilibration process of the quark-gluon plasma which have been undertaken 
over last two decades. In Sec. lllll various plasma instabilities are considered and the magnetic 
Weibel modes are argued to be relevant for the quark-gluon plasma produced in relativistic 
heavy-ion collisions. In Sec. IIVI I discuss how the unstable modes are initiated while in 



^ The term 'parton' is used to denote a quasiparticle fermionic (quark) or bosonic (giuon) excitation of the 
quark-gluon plasma. 

^ Although anisotropic systems are considered, the characteristic momentum in all directions is assumed to 
be of the same order. 
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Sec. Hthe mechanism of unstable mode growth is explained in terms of elementary physics. 
Sec. IVIl is devoted to solutions of the dispersion equation which provide dispersion relations 
of the unstable modes. In Sec. lVIIl it is explained why the instabilities isotropize the system. 
A phenomenon of spontaneous abelianization of the system's configuration is considered in 
the same section. The two next sections contain more formal material. The Hard Loop 
effective action of anisotropic plasma is presented in Sec. IVlllI while Sec. IIXI deals with 
the equations of motion which are used to study temporal evolution of anisotropic plasma. 
Results of recent numerical simulations of the plasma evolution are presented in Sec. |3 The 
review is closed with a brief discussion on possible signals of the instabilities and on desired 
improvements of theoretical approaches to the unstable quark-gluon plasma. 

Throughout the article there are used the natural units with h = c = ks = 1] the metric 
convention is (1, — 1, — 1, — 1); the coupling constant = g'^/4:H is assumed to be small; 
quarks and gluons are massless. 



II. EQUILIBRATION OF THE QUARK-GLUON PLASMA 



To present the scenario of instabilities driven isotropization in a broader context, I start 
with a brief review of numerous attempts to understand the equilibration processes of the 
quark-gluon. The problem was posed over twenty years ago when the real prospects to create 
the quark-gluon plasma in terrestrial experiments appeared. Already in the early papers 
pubhshed in the eig hties H H El El S 

l2ll . I22I I , main directions of further studies were 
drawn. The space-time structure of ultrarelativistic heavy-ion collisions was found [l9| to 
provide an estimate of the system's temperature and the lower bound of the thermalization 
time. The Boltzmann equation in the relaxation time approximation and the Fokker- 
Planck equation were used to follow the equilibration process. The Schwinger mechanism 
of particle production was included in kinetic theory treatment of the thermalization 3, 2^ 
and the pure perturbative mechanism was analysed as well '21^. The equilibration was also 
studied within the Monte Carlo parton cascade model which, however, took into account 
only binary parton-parton collisions. 

These lines of research were continued in the next decade. The parton cascade approach 
was greatly improved ^[ by, in particular, including the gluon radiation in the initial 
and final states of parton-parton interactions. The radiation proved to be very important 
for the equilibration process 0, . These detailed numerical studies are summarized in 
the review [2^. Another perturbative parton cascade approach combined with the string 
phenomenology for non-perturbative interactions is presented in |27|. The analytical studies 
of the thermalization were continued in 2^ 2^ 3^, 31, 3^, see also j27l|. where, in particular, 
the gluons were convincingly shown to equilibrate much faster than the quarks, the free 
streaming and the role of infrared cut-offs in the parton-parton cross sections were elucidated. 
Much efforts were invested in the studies of multi-particle processes js^ 0, which 



were already implemented in the parton cascade type models j26|, |27[ . The inelastic process 
2 3 attracted a lot of attention. Although it is of higher order in as, it is responsible for 
the parton number equilibration and it dominates the entropy production 3^ 3^, 3^ . 

There are two very recent transport theory approaches to the equilibration problem based 
on big numerical codes where the role of the multi-particle processes is emphasized (37i 
3^]. The authors of jl^] include particle production and absorption via the process 2 3 
while the three-particle collisions 3 3 are studied in |38,J. Within both approaches the 
equilibration is claimed to be significantly speeded-up when compared to the equilibration 
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driven by the binary collisions. However, the interaction rates of multi-particle processes 
are known to suffer from severe divergences, and thus, the actual role of the multi-particle 
interactions crucially depends on how the rates are defined, computed and regularized. 

The observation that the multi-particle interaction rates are sometimes divergent was 
actually used to explain the very fast equihbration of the quark-gluon system. The so-called 
collinear divergences of the gluon multiplication process 2^3 cancel in the equilibrium. If 
the cancellation does not occur in the non-equilibrium systems, as argued in |31|], the equi- 
libration, which is driven by very large - formally divergent - interaction rates, is extremely 
fast even in the weakly coupled plasma 139^ . 

The thermalization of the quark-gluon plasma was also discussed from a very different 
point of view where the equilibration is not due to the inter-parton collisions but due to the 
chaotic dynamics of the non-Abelian classical fields (coupled or not to the classical coloured 
particles) |4ol 41 1, see also a very recent paper i^]. Then, the equilibration time is controlled 



by the maximal Lyapunov exponent. 

At the turn of the millennium, when a large volume of experimental data from the 
RHIC started to flow, understanding of the equilibration process became a burning issue 
as the data favoured a very short equilibration mentioned in the Introduction. Within 
the concept of strongly coupled quark-gluon plasma, the problem is trivially solved as the 
strongly interacting system is indeed equilibrated very fast. However, it is still an open issue 
whether the plasma at the collision early stage is indeed strongly coupled. 

A novel development concerned a treatment of the initial state of the parton system 
which evolves towards equilibrium. In the papers mentioned above, one usually assumed 
that the initial partons are produced due to the (semi-)hard interactions of partons of the 
incident nuclei. Thus, jets and minijets form such an initial state which can be parametrized 
in several ways [i^. Recent studies of the equilibration problem which adopt the minijet 
initial conditions are presented in 3 S i3- 

In the already mentioned 'bottom- up' thermalization scenario the initial state was 
assumed to be shaped by the QCD saturation mechanism. Then, the initial state is dom- 
inated by the small x gluons of transverse momentum of order Qs which is the saturation 
scale. These gluons are freed from the incoming nuclei after a time Qj^- Weak coupling 
techniques are applicable as Qs is expected to be much smaller than Aqcd at sufficiently high 
collision energies. The saturation mechanism is incorporated in the effective field approach 
known as the Colour Glass Condensate ji^ where the small x partons of large occupation 
numbers are treated as classical Yang-Mills fields. Hard modes of the classical fields play 
the role of particles here. The equilibration processes with the minijet and saturation initial 
states were compared to each other in ji^ . 

The 'bottom-up' thermalization scenario where not only binary collisions but the 
processes 2 <-> 3 are included, takes into account the system's expansion. The equilibration 
processes splits into several stages parametrically characterized by a^Q~^ where n is a 
fractional power. The thermalization time is of order a~^^^^^Q^^ . However, as stressed 
in the Introduction, the collisional isotropization is apparently too slow to comply with 
the experimental data. The calculations performed within the 'bottom-up' scenario i(y\ were 
criticized pL^J for treating the parton momentum distribution as isotropic, and thus, ignoring 
the instabilities which actually speed up the equilibration process. Recently, an influence of 
the instabilities on the 'bottom- up' time scales has been discussed in ji^. It has been also 
argued js^] that a somewhat modified scenario remains valid for a sufficiently late stage of 
the equilibration process when the instabilities are no longer operative. 
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At the end I mention rather unconventional approaches to the fast equilibration problem. 
It was argued in 51, 5^ 5^ that the momentum distribution of partons is of the equilibrium 
form just after the production process. Thus, the very process of particle production leads 



to the equilibrium state without any secondary interactions. The authors of [51|, |52| refer 



to the Schwinger mechanism of particle's production due to the strong chromoelectric field. 
The transverse momentum but not longitudinal one is claimed to be 'equilibrated' in this 



way 



51,^ 



The key ingredient of the approach |53[, where the longitudinal momentum 
is also thermal, is the Hawking- Unruh effect: an observer moving with an acceleration a 
experiences the influence of a thermal bath with an effective temperature a/27r, similar to 
the one present in the vicinity of a black hole horizon. The idea behind the approaches 



51, B 



53l | is elegant and universal - it can be applied not only to nucleus-nucleus but to 
hadron-hadron or even to — e~ collisions - but it cannot explain how the equilibrium 



— e 

state is maintained when the parton's free streaming drives the system out of equilibrium. 
Secondary interactions are then certainly needed. 

Finally, I note a very interesting 'no-go' theorem [s^, 5^, which states that the pertur- 
bative thermalization is impossible, as any Feynman diagram of any order leads in the long 
time limit to the time scaling of the energy density corresponding to the free streaming not 
to the Bjorken hydrodynamics. However, it is not quite clear whether the theorem applies 
to the relativistic heavy-ion collisions as the equilibrium state of matter produced in the 
collisions is presumably only a transient state which changes into free streaming at the late 
times of the system's evolution. 



III. RELEVANT PLASMA INSTABILITIES 

The electron-ion plasma is known to experience a large variety of instabilities Those 
caused by coordinate space inhomogeneities, in particular by the system's boundaries, are 
usually called hydrodynamic instabilities, while those due to non-equilibrium momentum 
distribution of plasma particles are called kinetic instabilities. Hardly anything is known 
about hydrodynamic instabilities of the quark-gluon plasma, and I will not speculate about 
their possible role in the system's dynamics. The kinetic instabilities are initiated either 
by the charge or current fluctuations. In the first case, the electric field (E) is longitudinal 
(E II k, where k is the wave vector), while in the second case the field is transverse (E _L k). 
For this reason, the kinetic instabilities caused by the charge fluctuations are usually called 
longitudinal while those caused by the current fluctuations are called transverse. Since 
the electric field plays a crucial role in the longitudinal mode generation, the longitudinal 
instabilities are also called electric while the transverse ones are called magnetic. In the 
non-relativistic plasma the electric instabilities are usually much more important than the 
magnetic ones, as the magnetic effects are suppressed by the factor f^/c^ where v is the 
particle's velocity. In the relativistic plasma both types of instabilities are of similar strength. 
The electric instabilities occur when the momentum distribution of plasma particles has more 
than one maximum, as in the two-stream system. A sufficient condition for the magnetic 
instabilities is, as discussed in Sec. IVH anisotropy of the momentum distribution. 

Soon after the concept of quark-gluon plasma had been established, the existence of the 
colour kinetic instabilities, fuUv analogous to those known in the electro dynamic plasma, was 
suggested 57, IH, |5^ Ig^, 61, Ig^, 6^- In these early papers, however, there was considered 



a two-stream system, or more generally, a momentum distribution with more than one 
maximum. While such a distribution is common in the electron-ion plasma, it is rather 
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irrelevant for the quark-gluon plasma produced in relativistic heavy-ion collisions where the 
global as well as local momentum distribution is expected to monotonously decrease in every 
direction from the maximum. The electric instabilities are absent in such a system but, as 
demonstrated in {64] , a magnetic unstable mode known as the filamentation or Weibel 
instability jH^ is possible. The filamentaion instability was shown [g^I to be relevant 
for the quark-gluon plasma produced in relativistic heavy-ion collisions as the characteristic 
time of instability growth is shorter or at least comparable to other time scales of the parton 
system evolution. And the instabilities - usually not one but several modes are generated 
- drive the system towards isotropy, thus speeding up its equilibration. In the following 
sections a whole scenario of the instabilities driven equilibration is reviewed. 

IV. SEEDS OF FILAMENTATION 

Let me start with a few remarks on degrees of freedom of the quark-gluon plasma. Various 
problems will be repeatedly discussed in terms of classical fields and particles which are only 
approximate notions in the quark-gluon plasma being a system of relativistic quantum fields. 
However, collective excitations, which are bosonic and highly populated, can be treated as 
classical fields while bosonic or fermionic excitations, with the energy determined by the 
excitation momentum (due to the dispersion relation), can be treated as (quasi-)particles. 
In the weakly coupled quark-gluon plasma in equilibrium, an excitation is called hard when 
its momentum is of order T, which is the system's temperature, and it is called soft when 
its momentum is of order gT. Within the Hard Loop dynamics, the hard excitations can be 
treated as particles while the gluonic soft excitations as classical fields [(36]. It is expected 
that a similar treatment is possible in the non-equilibrium plasma as well. Thus, the terms 
partons, quarks, gluons, particles will be used to denote quasiparticle hard excitations. The 
classical chromodynamic field will represent gluonic soft collective excitations. 

After the introductory remarks, let me discuss how the unstable transverse modes are 
initiated. For this purpose I consider a parton system which is homogeneous but the parton 
momentum distribution is not of the equilibrium form, it is not isotropic. The system is on 
average locally colourless but colour fluctuations are possible. Therefore, {ja{x)) = where 
ja{x) is a local colour four-current in the adjoint representation of SU(A^c) gauge group with 
/i = 0, 1, 2, 3 and a = 1,2, ... , {N^ — 1) being the Lorentz and colour index, respectively; 
X = (t, x) denotes a four-position in coordinate space. 

Since I assume that the quark-gluon plasma is weakly coupled, the non-interacting gas of 
quarks, antiquarks and gluons can be treated as a first approximation. As discussed in detail 
in (g^, the current correlator for a classical system of non-interacting massless partons is 

Mr(t,x) {ji^{t„^,)jah,^,)) = '^"7 ^ ^^p^ 5(^)(x- vt) , (1) 

where v = p/|p|, (t,x.) = (^2 — ti,^2 — ^i) and the effective parton distribution function 
/(p) equals n{p) + n(p) + 2Ncng{p); n{p), n(p) and ng{p) give the average colourless 
distribution function of quarks Q^^{p,x) = 6^^n{p), antiquarks Q^^{p,x) = (5*%(p), and 
gluons G°'^{p,x) = 6'^^ng{p). The distribution function of (anti-)quarks and gluons are 
matrices belonging to the fundamental and ajoint representation, respectively, of the SU(A^'c) 
gauge group. Therefore, i, j = 1,2, Nc and a,b = 1, 2, {N^ — 1). 

Due to the average space-time homogeneity, the correlation tensor (^J depends only on 
the difference (^2 — ii,X2 — xi). The space-time points (ti,xi) and (^25X2) are correlated 
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in the system of non-interacting particles if a particle travels from (ti,xi) to (t2,X2). For 
this reason the delta (5*-^^(x — vt) is present in the formula (^. The momentum integral of 
the distribution function simply represents the summation over particles. The fluctuation 
spectrum is found as a Fourier transform of the tensor i.e. 

Mr(^, k) = ^ / '^'^7 (Cf^ ^ ^^p^ 2^'^^^ " ^""^ ■ 

To compute the fluctuation spectrum, the parton momentum distribution has to be spec- 
ified. Such calculations with two forms of the momentum distribution are presented in 
[gtI ]. Here I only qualitatively discuss Eqs. ((TJ |2)), assuming that the parton momentum 
distribution is anisotropic. 

In heavy-ion collisions, the anisotropy is a generic feature of the parton momentum dis- 
tribution in a local rest frame. After the first collisions, when the partons are released from 
the incoming nucleons, the momentum distribution is strongly elongated along the beam - 
it is of the prolate shape with the average transverse momentum being much smaller than 
the average longitudinal one. Due to the free streaming, it evolves in the local rest frame 
to the distribution which is squeezed along the beam - it is of the oblate shape with the 
average transverse momentum being much larger than the average longitudinal one. In most 
cases, I assume that the distribution is elongated along the z axis but my considerations 
remain valid for the distribution, which is squeezed along the z axis, but the axes should be 
relabeled. 

With the momentum distribution elongated in the z direction, Eqs. (^Q ^ clearly show 
that the correlator M^^ is larger than M^^ or M^^. It is also clear that M^^ is the largest 
when the wave vector k is along the direction of the momentum deficit. Then, the delta 
function 6{uj — kv) does not much constrain the integral in Eq. Since the momentum 
distribution is elongated in the z direction, the current fluctuations are the largest when the 
wave vector k is the x — y plane. Thus, I conclude that some fluctuations in the anisotropic 
system are large, much larger than in the isotropic one. An anisotropic system has a natural 
tendency to split into the current filaments parallel to the direction of the momentum surplus. 
These currents are seeds of the filamentation instability. 



V. MECHANISM OF FILAMENTATION 

Let me now explain in terms of elementary physics why the fluctuating currents, which 
flow in the direction of the momentum surplus, can grow in time. To simplify the discus- 
sion, which follows [67], I consider an electromagnetic anisotropic system. The form of the 
fluctuating current is chosen to be 

j(a;) = j cos(fc^x) , (3) 

where is the unit vector in the z direction. As seen in Eq. there are current filaments 
of the thickness 7r/\kx\ with the current flowing in the opposite directions in the neighbouring 
filaments. 

The magnetic field generated by the current Q is given as 

= j-Gy sm{k^x) , 
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FIG. 1: The mechanism of filamentation insta- 
bihty, see text for a description. 
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and the Lorentz force acting on the partons, which fly along the z direction, equals 

j 

F{x) = g V X B(x) = -qv^ — e^. sin(A;^x) , 

where q is the electric charge. One observes, see Fig. [TJ that the force distributes the partons 
in such a way that those, which positively contribute to the current in a given filament, are 
focused in the filament centre while those, which negatively contribute, are moved to the 
neighbouring one. Thus, the initial current is growing and the magnetic field generated by 
this current is growing as well. The instability is driven by the the energy transferred from 
the particles to fields. More specifically, the kinetic energy related to a motion along the 
direction of the momentum surplus is used to generate the magnetic field. The mechanism 
of Weibel instability is explained somewhat differently in . 



VI. DISPERSION EQUATION 



The Fourier transformed chromodynamic field A'^{k) satisfies the equation of motion as 

'eg^"" -k^'k" -n>'''{k)]A,{k) = , (4) 

where k = {uj, k) and Il'^'^{k) is the polarization tensor or gluon self-energy which is discussed 
later on. Since the tensor is proportional to a unit matrix in the colour space, the colour 
indices are dropped here. A general plasmon dispersion equation is of the form 



det 



. 



(5) 



Equivalently, the dispersion relations are given by the positions of poles of the effective 
gluon propagator. Due to the transversality of Il'^^{k) {k^Il'^^{k) = k^Il^'^(k) = 0) not all 
components of Il^'^{k) are independent from each other, and consequently the dispersion 
equation which involves a determinant of a 4 x 4 matrix, can be simplified to the 
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determinant of a 3 x 3 matrix. For this purpose I introduce the colour permittivity tensor 
e^^ik) where the indices l,m,n = 1,2,3 label three-vector and tensor components. Because 
of the relation 

e'^{k)E\k)E^{k) = U^''{k)A^{k)A,{k) , 

where E is the chromoelectric vector, the permittivity can be expressed through the polar- 
ization tensor as 



Then, the dispersion equation gets the form 



det 



2Jm/ 







(6) 



The relationship between Eq. ^ and Eq. (jH)) is most easily seen in the Coulomb gauge when 
= and k • A(A;) = 0. Then, E = iuA and Eq. (j3]) is immediately transformed into an 

equation of motion of E(A;) which further provides the dispersion equation (jHj). 

The dynamical information is contained in the polarization tensor Il'^'^{k) given by 

Eq. (Illj) or, equivalently, in the permittivity tensor e'™(/c) which can be derived either 



within the transport theory or diagrammatically j68[. The result is 



e™(cu,k) = (5" 



+ — 

2u 



(27r)3a;-kv + i0+ dp^ 



l„,m 



kv\ fc'f 

1 U'™ H 



UJ 



(7) 



As already mentioned, the colour indices are suppressed here. 

Substituting the permittivity ((Zj) into Eq. (jHl), one fully specifies the dispersion equation 
(jni) which provides a spectrum of quasi-particle bosonic excitations. A solution co'(k) of 
Eq. (jni) is called stable when Imcj < and unstable when \m.uj > 0. In the first case the 
amplitude is constant or it exponentially decreases in time while in the second one there is 
an exponential growth of the amplitude. In practice, it appears difficult to find solutions 
of Eq. (jni) because of the rather complicated structure of the tensor (|7j). However, the 
problem simplifies as we are interested in specific modes which are expected to be unstable. 
Namely, we look for solutions corresponding to the fiuctuating current in the direction of 
the momentum surplus and the wave vector perpendicular to it. 

As previously, the momentum distribution is assumed to be elongated in the z direction, 
and consequently the fluctuating current also flows in this direction. The magnetic fleld has 
a non-vanishing component along the y direction and the electric fleld in the z direction. 
Finally, the wave vector is parallel to the axis x, see Fig. It is also assumed that the 
momentum distribution obeys the mirror symmetry /(— p) = /(p), and then the permit- 
tivity tensor has only non-vanishing diagonal components. Taking into account all these 
conditions, one simplifles the dispersion equation ® to the form 



H(uj) 



(u;, k^) = 



(8) 



where only one diagonal component of the dielectric tensor enters. 

It appears that the existence of unstable solutions of Eq. (jH)) can be proved without 
solving it. The so-called Penrose criterion j|56|], which follows from analytic properties of the 
permittivity as a function of u, states that the dispersion equation H{uj) = has unstable 
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solutions if H{uj = 0) < 0. The Penrose criterion was apphed to the equation JHl in 
but a much more general discussion of the instability condition is presented in [1311. Not 
entering into details, there exist unstable modes if the momentum distribution averaged 
(with a proper weight) over momentum length is anisotropic. 

To solve the dispersion equation (jS)), the parton momentum distribution has to be speci- 
fied. Several analytic (usually approximate) solution of the dispersion equation with various 
momentum distributions can be found in 0, 3, 3, 6^ • A typical example of the numerical 



solution, which gives the unstable mode frequency in the full range of wave vectors is shown 



in Fig. 121 taken from ll| . The momentum distribution is of the form 



^^P^^ (p2 +ai)3^ " ' 

where p± = ^pl + Py. The mode is pure imaginary and = liauj(k±). The value of the 

coupling is = fi'^/47r = 0.3, a± = 0.3 GeV and the effective parton density is chosen to 
be 6 fm"l As seen, there is a finite interval of wave vectors for which the unstable modes 
exist. 

The dispersion equation (jH} corresponds to a simple configuration where the wave vector 
is parallel to the axis x and it points to the direction of the momentum deficit while the 
chromoelectric field is parallel to the axis z and it points to the momentum surplus. However, 
there are more general unstable modes which are not aligned along the symmetry axes of 
the momentum distribution of particles. The wave vectors k and chromoelectric fields E 
of these modes have non-vanishing components in the directions of the momentum deficit 
and momentum surplus, respectively, and E is no longer perpendicular to k. Such unstable 
modes are discussed in llll. A quite general analysis of the dispersion equation of anisotropic 
systems is given in |l2l 169^. There is considered a class of momentum distributions which 
can be expressed as 

/(p) = /iso(Vp^ + e(np)2) , (9) 

where /iso(|p|) is an arbitrary (isotropic) distribution, the unit vector n defines a preferred 
direction and the parameter C, € (—1, c>o) controls the magnitude of anisotropy. 

As explained above, the existence of the unstable gluonic modes is a generic feature of 
the anisotropic plasma - even a weak anisotropy generates the instability. In contrary, the 
quark modes seem to be always stable. Although, a general proof of the quark mode stability 
is lacking, the modes appear to be stable even in the case of extremely anisotropic parton 
momentum distribution as in the two-stream system [l^l • Presumably, the quark modes are 
always stable because their population is constrained by Pauli blocking ^7^. 



VII. ISOTROPIZATION AND ABELIANIZATION 

When the instabilities grow the system becomes more isotropic because the Lorentz force 
changes the particle's momenta and the growing fields carry an extra momentum. To explain 
the mechanism I assume, as previously, that initially there is a momentum surplus in the 
z direction. The fluctuating current flows in the z direction with the wave vector pointing 
in the x direction. Since the magnetic field has a y component, the Lorentz force, which 
acts on partons flying along the z axis, pushes the partons in the x direction where there 
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FIG. 3: The effective potential of the unstable magnetic mode as a function of magnitude of two 
colour components of A*^ belonging to the SU(2) gauge group. The figure is taken from [t^. 



is a momentum deficit. Numerical simulations discussed in Sec. IXIshow that growth of the 
instabilities is indeed accompanied with the system's fast isotropization. 

The system isotropizes not only due to the effect of the Lorentz force but also due to the 
momentum carried by the growing field. When the magnetic and electric fields are oriented 
along the y and z axes, respectively, the Poynting vector points in the direction x that is 
along the wave vector. Thus, the momentum carried by the fields is oriented in the direction 
of the momentum deficit of particles. 

Unstable modes cannot grow to infinity and even in the electron-ion plasma there are 
several possible mechanisms which stop the instability growth [t^]. The actual mechanism 
depends on the plasma state as well as on the external conditions. In the case of the quark- 
gluon plasma one suspects that non-Abelian non-linearities can play an important role here. 
An elegant argument [t^ suggests that the non-linearities do not stabilize the unstable 
modes because the system spontaneously chooses an Abelian configuration in the course of 
instability development. Let me explain the idea. 

In the Coulomb gauge the effective potential of the unstable configuration has the form 

which is shown in Fig. El taken from [l^ . The first term (with > 0) is responsible 
for the very existence of the instability. The second term, which comes from the Yang- 
Mills lagrangian, is of pure non-Abelian nature. The term appears to be positive and thus 
it counteracts the instability growth. However, the non-Abelian term vanishes when the 
potential A" is effectively Abelian, and consequently, such a configuration corresponds to 
the steepest decrease of the effective potential. Thus, the system spontaneously abelianizes 
in the course of instability growth. In Sec. [Xj where the results of numerical simulations are 
presented, the abelianization is further discussed. 

VIII. HARD-LOOP EFFECTIVE ACTION 

Knowledge of the gluon polarization tensor or, equivalently, the chromoelectric permit- 
tivity tensor is sufficient to discuss the system's stability and the dispersion relations of 
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unstable modes. For more detailed dynamical studies the effective action of anisotropic 
quark-gluon plasma is needed. Such an action for a system, which is on average locally 
colour neutral, stationary and homogeneous, was derived in The starting point was 

the effective action which describes an interaction of classical fields with currents induced 
by these fields in the plasma. The lagrangian density is quadratic in the gluon and quark 
fields and it equals 

C,{x) = -J d'y(^^A';{x)U^,a^ - y)Kiy) + - y)^{y)^ ; (10) 

the Fourier transformed gluon polarization tensor 11^^ (A;) and the quark self-energy S(A;) 
read 

-2 " d^p /(p) {p ■ k){kf'p'' + pf^k") - k'^p^'p" - {p ■ kfg'"' 



Oab — 



2 J (27r)3 IpI 



(p ■ kf 



8Nr 



d^p /(p)p-7 

{2TTy IpI p ■ k 



(12) 



where /(p) and /(p) are the effective parton distribution functions defined as /(p) = n{p) + 
n(p) + 2Ncng{p) and /(p) = n(p) + n(p) + 2ng(p); n(p), n(p) and ng{p) are, as already 
mentioned below Eq. ((T)), the distribution functions of quarks, antiquarks and gluons of single 
colour component in a homogeneous and stationary plasma which is locally and globally 
colourless; the spin and flavour are treated as parton internal degrees of freedom. The 
quarks and gluons are assumed to be massless. The polarization tensor (jTTll can be derived 
within the semiclassical transport theory 0, |^ or diagrammatically [68], following the 
formal rules of the Hard Thermal Loop app roach. The quark self-energy (fT^ has been 
derived so far only diagrammatically 68 but the derivation is also possible within the 
transport theory, as it has been done in 76^ for the case of equilibrium plasma. The action 
(fTUj) holds under the assumption that the field amplitude is much smaller than T/g where 
T denotes the characteristic momentum of (hard) partons. 

Following Braaten and Pisarski Q| , the lagrangian pO|) was modified to comply with the 
requirement of g auge invariance. The final result, which is non-local but manifestly gauge 
invariant, is UA 



9^ f 



d^p 
(2^ 



/(P) 



/ pUpP 



\{p-D) 



p ' 



x) 



(13) 



+ i- 



P-1 



ANr 



/(p) m{x)^^>{x 
p ■ D 



where F^*^ is the strength tensor and D denotes the covariant derivative. The effective 
action ()13p generates n— point functions which obey the Ward-Takahashi identities. For 
the equilibrium plasma the action (1131) is equivalent to that one derived in [78^1 and in the 
explicitly gauge invariant form in [77[. The equilibrium Hard Loop action was also found 
within the semiclassical kinetic theory 



IX. EQUATIONS OF MOTION 



Transport theory provides a natural framework to study temporal evolution of non- 
equilibrium systems and it has been apphed to the quark-gluon plasma for a long time. 



12 



The distribution functions of quarks (Q), antiquarks (Q), and gluons (G), which are the 



the form SOUSlI: 



Nc X Nc and (N? — 1) x (A^^ — 1) matrices, respectively, satisfy the transport equations of 

dQ{p,x) 



dpu 
dQ{p,x) 

dpu 
dG{p, x) 

dpu 



0, 











(14) 



where {...,...} denotes the anticommutator; the transport equation of (anti-) quarks is written 
down in the fundamental representation while that of gluons in the adjoint one. Since the 
instabilities of interest are very fast, much faster than the inter-parton collisions, the collision 
terms are neglected in Eqs. (fT^ . The gauge field, which enters the transport equations (dlj), 
is generated self-consistently by the quarks and gluons. Thus, the transport equations (fT^ 
should be supplemented by the Yang-Mills equation 



D.F^'ix) 



fix) 



where the colour current is given as 



J^ix) = -g 



p^ 



d^p 



Tr[r,(g(p,x) -g(p,x))] +Tr[T,G(p,x)] 



(15) 



(16) 



with Ta and being the SU(Ai'c) group generators in the fundamental and adjoint repre- 
sentation, respectively. There is a version of the equations fjl4t [T^ where colour charges of 
partons are treated as a classical variable [82|. Then, the distribution functions depend not 
only on x and p but on the colour variable as well. 

When the equations (|141 ITH|) are linearized around the state, which is stationary, homoge- 
neous and locally colourless, the equations provide the Hard Loop dynamics encoded in the 
effective action (fT^ . The equations are of particularly simple and elegant form when the 
quark SQ{p, x), antiquark 6Q{p, x) and gluon 6G{p, x) deviations from the stationary state 
described by Qo (p) = 5*%(p), Qo (p) = '^*''^(p)) Gq^{p) = 6°''^ng{p) are parameterised 
by the field iy^(v,x) through the relations 



6Q{p,x) 



g^W^i.,x) 



dpt^ 



6G{p,x) = g 



5Qip,x) = -g^^W>^{v,x) 



dpi' 



dpf^ 



r,Tr[r,l^^(v,x)] 



where v = p/|p|. Then, instead of the three transport equations (fTH) one has one equation 



v^D>'W''{^r,x) 
while the Yang-Mills equation ()15|) reads 



D^F^^'ix) 



-g 



-VpFP'^ix) 



<fp p^ df{p) 



(27r)3 IpI dpp 



^^(v,^) 



(17) 



(18) 
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FIG. 4: Time evolution of the (scaled) energy den- 
sity (split into various electric and magnetic com- 
ponents) which is carried by the chromodynamic 
field. The simulation is 1+1 dimensional and the 
gauge group is SU(2). The parton momentum dis- 
tribution is squeezed along the z axis. The solid 
line corresponds to the total energy transferred 
from the particles. The figure is taken from [l3 |. 
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FIG. 5: Time evolution of the kinetic 
energy of particles (upper panel) and of 
the energy of electric and magnetic fields 
(lower panel) in GeV/fm'^ for the U(l) 
and SU(2) gauge groups. The figure is 
taken from llSll. 



where = (1, v) and, as previously, /(p) = n{p) + ^(p) + 2Ncng{p). In contrast to the 
effective action (fT^ . the equations (fTTl ITS]) are local in coordinate space. Therefore, the 
transport equation ()17p combined with Eq. ()18|) is often called local representation of the 
Hard Loop dynamics. The equations ()171 ITHj) . which for the isotropic eq uilibrium plasma 
were first given in [g^, are used in the numerical simulations discussed in the 

next section. 



X. NUMERICAL SIMULATIONS 



Temporal evolution of the aniso trop ic quark-gluon plasma has been recently studied by 
means of numerical simulations 3 3 lEIi^- The simulations, which have been performed 
in two very different dynamical schemes by three groups of authors, are of crucial importance 
as they convincingly demonstrate a key role of the instabilities in the evolution of anisotropic 
quark-gluon plasma. 

The dynamics governed by the Hard Loop action (fT^ and described by the equations 
(fTTl ITS|) has been simulated in l^, 8^ "EM- These simulations provide fully a reliable infor- 
mation on the field dynamics provided the potential's amplitude is not too large: A'^ <^ T/g 
where T is the characteristic momentum of (hard) partons. Since the equations (fTTl ITHj) de- 
scribe small deviations from the stationary homogeneous state, only a small fraction of the 
particles is influenced by the growing chromodynamic fields. Therefore, the (hard) particles 
effectively play a role of the stationary (anisotropic) background. In the simulation the 
classical version of the equations (jl41 IT^ is used. The quark-gluon plasma is treated as a 
completely classical system: partons, which carry classical colour charges, interact with the 
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FIG. 6: Temporal evolution of the functionals 
C and (^i-ins measured in GeV. The figure is 
taken from llSll. 
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FIG. 7: Temporal evolution of the (scaled) 
functionals C and jrms- The figure is taken 
from (I4l |. 



self-consistently generated classical chromo dynamic field. 

The simulations 3 3] have been effectively performed in 1+1 dimensions as the chro 



modynamic potentials depend on time and one space variable. The calculations 
represent full 1+3 dimensional dynamics. In most cases the SU(2) gauge group was studied 
but some SU(3) results, which are qualitatively very similar to 31/(2) ones, are given in (s^ . 

The techniques of discretization used in the simulations j3, EE EE 13 are rather differ- 
ent while the initial conditions are quite similar. The initial field amplitudes are distributed 
according to the Gaussian white noise and the momentum distribution of (hard) partons is 
strongly anisotropic. For example, in the classical simulation [l^ the initial parton momen- 
tum distribution is chosen as 



/(p) ~ KVx) e f'h-d , (19) 

with phard = 10 GeV. The results are actually insensitive to the specific form of the momen- 
tum distribution. If the parton distribution function is written in the form Q, the results 
are shown 0, ^ to depend only on two parameters: ^ and the Debye mass rriD of the 



corresponding isotropic system i.e. 



47r^ JO dp 

In Fig. taken from 3, the results of the Hard-Loop simulation performed in 1+1 



dimensions are shown. One observes exponential growth of the energy density stored in 
fields and the energy density is dominated, as expected, by the magnetic field which is 
transverse to the direction of the momentum deficit. The growth rate of the energy density 
appears to be equal to the growth rate 7* of the fastest unstable mode. Fig. El which is taken 
from [1^, shows results of the classical simulation on the 1+1 dimensional lattice of physical 
size L = 40 fm. As in Fig. EJ the amount of field energy, which is initially much smaller 
than the kinetic energy of all particles, grows exponentially and the magnetic contribution 
dominates. 

The Abelian (U(l)) and non-Abelian (SU(2)) results of the 1+1 dimensional simulation 
presented in Fig. El are remarkably similar to each other. The abelianization, explained in 
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FIG. 8: Time evolution of the (scaled) 
chromomagnetic energy density in the 
1+3 dimensional simulation. The 
Abelian result and that of 1+1 dimen- 
sions are also shown. The figure is taken 
from l83l. 




40 50 



FIG. 9: Time evolution of the (scaled) energy den- 
sity (split into various electric and magnetic com- 
ponents) of the chromodynamic field in the 1+1 
and 1+3 simulations. 'HL' denotes the total en- 
ergy contributed by hard particles. The figure is 
taken from kJI. 



Sec. IVm appears to be very efficient in 1+1 dimensions, as shown in Figs. El [71 taken from 



isj l and 1J|, respectively. The authors of ^] analysed the functional 



/ rL dx - rfx\/Tr[(iK,AJ)21 



which were introduced in [73[ . The quantities jrms and C, studied in jlj| and shown in Fig. [71 
are fully analogous to 0rms and C defined by Eq. (^0]) but the components of chromodynamic 
potential are replaced by the respective components of colour current. As seen in Figs. IHIIZI 
the field (current) commutator decreases in time although the magnitude of field (current), 
as quantified by 0rms (jrms), grows. 

It is worth mentioning that the functionals (|2m) defined through the gauge potentials are 
gauge invariant provided the potentials depend only of one time and one space variables and 
the gauge transformations preserve this property. Thus, the functionals are well suited 
for 1 + 1 dimensional simulations. However, the functionals (|2m) are not gauge invariant 
under general 1 + 3 dimensional gauge transformations. When the potential components 
are replaced by the respective current components, as proposed in 1J|, the functionals are 



gauge invariant not only under 1 + 1 but also under 1 + 3 dimensional transformations. 

The results of the 1+3 dimensional simulations f83','83] appear to be qualitatively different 
from those of 1+1 dimensions. As seen in Figs. IHlini taken from and ^], respectively, 
the growth of the field energy density is exponential only for some time, and then the 
growth becomes approximately linear. It appears that the regime changes when the field's 
amplitude is of order k/g where k is the characteristic field wave vector. Then, the non- 
Abelian effects start to be important. Indeed, Fig.^1 which is taken from ^8^, demonstrates 
that the abelianization is efficient in the 1+3 dimensional simulations [83, Mj only for a finite 
interval of time. The commutator C shown in Fig. ^1 which is a natural generalization of 
the 1+1 dimensional commutator defined by Eq. (j^Uj) with the current components instead 
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FIG. 10: Temporal evolution of the field 
commutator quantified by C. The figure is 
taken from 18311 . 



10=^ 



in 10 
i_ 

3 
CO 
CO 

£ 1 

Q. 
O 



10' 



0.2 



g T,,, SU(2) 

g (T,y+T,,)/2, SU(2) 

— gT,„u(i) 

g'(VTJ/2,U(1) 

o!4 ole ols 

time t/L/N, 



FIG. 11: Temporal evolution of the energy- 
momentum tensor components T^^ and {T^^ + 
r^^)/2. The Abelian and non-Abelian results 
are shown. The figure is taken from [l5| . 



of the potential ones, first decreases but after some time it starts to grow and returns to its 
initial value. 

As discussed in two very recent papers jH, E^, the physics of the late stage of instability 
development when the energy stored in the fields grows linearly with time is very interesting. 
Similarly to the turbulence, the unstable modes do not grow any more but due to non- 
Abelian interactions the energy provided by the particles is cascaded towards harder and 
harder modes. 

The effect of isotropization due to the action of the Lorentz force is nicely seen in the 
1+1 dimensional classical simulation 1^|. In Fig. ^2 which is taken from [3], there are 
shown diagonal components of the energy-momentum tensor 



rpfiu 



(27r)3 



/(P) 



The initial momentum distribution is given by Eq. (fT^ . and consequently T^^ = at 
t = 0. As seen in Fig. T^^ exponentially grows. However, a full isotropy, which requires 

rpxx ^ ^j,yy ^ j.,,^ achieved. 



XI. OUTLOOK AND FINAL REMARKS 

One wonders whether the presence of the instabilities at the early stage of relativistic 
heavy-ion collisions is experimentally observable. The accelerated equilibration is obviously 



very important though it is only an indirect signal. It has been suggested [87|,|88[ that strong 
chromomagnetic fields generated by the instabilities can lead to a specific pattern of jet's 
deflections. This promising proposal, however, requires further studies. 

Another idea has been formulated in [s^. The quark-gluon plasma, which is initially 
anisotropic, is isotropized fast due to the magnetic instabilities. Such a non-equilibrium 
plasma manifests, as recently observed an approximate hydrodynamic behaviour even 
before the equilibrium is reached. The point is that the structure of the ideal fluid energy- 
momentum tensor i.e. T^^ = {e + p) u^u^ — pg^'^, where e, p and u'^ is the energy density. 
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pressure and hydro dynamic velocity, respectively, holds for an arbitrary but isotropic mo- 
mentum distribution, e and p are then not the energy density and pressure but the moments 
of the distribution function which are equal the energy density and pressure in the equilib- 
rium limit. Since the tensor T^^ always obeys the continuity equation d^T^'^ = 0, one gets 
an analogue of the Euler equation. However, due to the lack of thermodynamic equilibrium 
there is no entropy conservation and the equation of state is missing. 

The azimuthal fluctuations have been argued 8^| to distinguish the approximate hydro- 
dynamics - characteristic for the instabilities driven isotropization - from the real hydro- 
dynamics describing a system which is in a local thermodynamic equilibrium, as advocated 
by proponents of the strongly coupled plasma f3|]. Non-equilibrium fluctuations are usually 
significantly larger than the equilibrium fluctuations of the same quantity. A specific exam- 
ple of such a situation is given in Sec. IIVI where the current fluctuations in the anisotropic 
system are discussed. Thus, one expects that the (computable) fluctuations of V2 produced 
in the course of real hydrodynamic evolution are significantly smaller than those generated 
in the non-equilibrium quark-gluon plasma which is merely isotropic. It should be stressed 
here that the elliptic flow is generated in the collision relatively early stage when there is 
a large configuration-space asymmetry of the colliding system. Since a measurement of v 2 
fluctuations is rather difficult, it was also argued that an integral measurement of the 
azimuthal fluctuations can help as well to distinguish the equilibrium from non-equilibrium 
fluctuations. 

Further suggestions of detectable signals of the instabilities are very much needed. How- 
ever, an experimental verification will certainly require much better theoretical understand- 
ing of the equilibration process. Although an impressive progress has been recently achieved, 
the numerical simulations 3 3 SJ| ^^^^ still quite far from the situation which occurs 
in relativistic heavy-ion collisions. A realistic simulation should properly take into account 
an initial state of colliding nuclei; the system's expansion, which slows down or even cuts 
off the instabilities growth [13 , needs to be included; the simulation should be 1 + 3 
dimensional as the dimensionality really matters (ssl . 83 | ; the effect of back reaction of fields 
on the particles has to be taken into account to observe the system's isotropization. 

Some of the above requirements are met by a very recent numerical study [90. . 91i] where 
the initial state corresponds to the Colour Glass Condensate ji^ where small x partons 
of large occupation numbers, which dominate the wave functions of incoming nuclei, are 
treated as classical Yang-Mills fields. As already mentioned, hard modes of the classical 
fields play the role of particles here. The study shows that the instabilities, identified as 
the Weibel modes, are indeed generated when the system of Yang-Mills fields representing 
colliding nuclei expands into the vacuum. However, the unstable mode growth gued 
11, IJ, slowed down. 
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Understanding of the late stage of the instability growth, when fields are of large mag- 
nitude, is a real theoretical challenge. The mechanism of instab ility saturation is not well 
known even in the electron-ion plasma, see e.g. a recent paper Non-linear effects, in 



particular those of non-Abelian nature, are then essential. Except for the classical simula- 
tions [3, 0, im , the evolution of anisotropic quark-gluon plasma has been studied within 
the Hard Loop approximation. An attempt to go beyond it has been undertaken in 0| 
where the higher order terms of the effective potential of the anisotropic system have been 
computed. Since these terms can be negative, the instability is then driven not only by 
the negative quadratic term but by the higher order terms as well. However, before a real 
progress in the strong field domain can be achieved, one still needs a better insight into the 
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Hard Loop dynamics which has appeared to be very rich HHQ. 

In summary, the magnetic instabihties provide a plausible explanation of the surprisingly 
short equilibration time observed in relativistic heavy-ion collisions. The explanation does 
not require a strong coupling of the quark-gluon plasma. Fast isotropization of the system is 
a distinctive feature of the instabilities driven equilibration. Two signals of the instabilities 
have been suggested but quantitative predictions are lacking. New ideas are certainly needed. 
In spite of the impressive progress, which has been achieved recently, a theoretical description 
of the unstable quark-gluon plasma requires further improvements. 
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